Soliton Confinement and the Excitation Spectrum of Spin-Peierls Antiferromagnets * 
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The excitation spectrum of spin-Peierls antiferromagnets is discussed taking into acount phonon 
dynamics but treating inter-chain elastic couplings in mean field theory. This gives a ladder of 
soliton-anti-soliton boundstates, with no soliton continuum, until soliton deconfinement takes place 
at a transition into a non-dimerized phase. 
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Much of the theoretical work on spin-Peierls systems largely ignores phonon dynamics, regarding the lattice dis- 
tortion as being static, thus producing an alternating Heisenberg exchange coupling. In the approach of Cross and 
Fisher Q the resulting susceptibility of the spin system (its response to a lattice distortion) is then fed into an RPA 
calculation of the phonon Green's function. However, the magnetic excitations, at least in the dimerized phase, are 
regarded as containing no phonon component. On the other hand, some theoretical work on the magetic field or im- 
purity induced undimerized phase has considered solitons which are at the same time magnetic excitations (they have 
spin 1/2) and involve lattice dynamics (the lattice distortion switches between the two phases at the location of the 
soliton). Khomskii et al. have developed a simple appealing picture of the spin-Peierls transition (as a function of 
field, impurity concentration or temperature) based on soliton unbinding. Their approach differs fundamentally from 
I ' that of Cross and Fisher in its treatment of phonon dynamics. Khomskii ct al. attempt to treat the one-dimensional 
, phonon dynamics accurately, while approximating the interchain elastic couplings by mean field theory. This ap- 
proach probably works best when the transverse phonon dispersion energy is small compared to the magnetic energy 
gap. Such an approach leads to an effective one-dimensional model containing both spontaneous dimerization due 
to one-dimensional phonons and explicit dimerization produced by the mean field of the neighbouring chains. This 
is to be contrasted with the other approach which only contains explicit dimerization. In the approach of Khomskii 
f^. ' et al. the fundamental excitations of the system are solitons. In the spin-Peierls ordered phase the solitons (s) and 
Q\ , anti-solitons (s) are bound together in pairs by a linear potential arising from the explicit dimerization potential of 
the neighbouring chains. A sufficient temperature, field or impurity concentration drives this self-consistently deter- 
mined dimerizing field to zero, eliminating the linear potential between solitons, allowing free solitons to propogate. 
From this perspective the spin-Peierls transition corresponds essentially to soliton deconfinement rather than phonon 
softening. 

In this paper we wish to extend this approach to a more quantum mechanical treatment of the excitations in the spin- 
Peierls ordered phase. The ss system, with its linear potential is quantized, leading to a ladder of boundstates which 
can have spin or 1. The s=l boundstates correspond to magnons. As the system becomes more one-dimensional, the 
linear potential gets weaker and the number of stable boundstates increases, diverging in the one-dimensional limit. 
The s=0 boundstates could also be interpreted as boundstates of two magnons, but in the highly one-dimensional 
case they are better interpreted as weakly bound singlet ss pairs with energies given by twice the energy of the soliton 
plus a (positive) "binding energy" associated with the linear potential. All boundstates are both magnetic and elastic 
in character; the s=0 boundstates are not neccessarily distinct from optical phonons. 

Within this approach, a soliton continuum of excitations does not occur in the ordered phase, contrary to the claims 
in ||. Instead the continuum is quantized into a ladder of boundstates. A magnon continuum can occur, beginning 
at precisely twice the magnon gap. All excited ss boundstates must lie below this continuum (in the appropriate spin 
channel) in order to be stable. The occurance of ss boundstates in Heisenberg models with competing spontaneous and 
explicit dimerization was mentioned by Haldane [Q. The fact that a soliton continuum cannot occur in spin-Peierls 
systems due to soliton confinement was pointed out recently by Uhrig and Schulz 

The occurence of a ladder of ss boundstates is generic to quasi-one-dimensional systems with broken discrete 
symmetries, at least within this type of mean field treatment. A very similar approach was taken by Shiba jf| to 
Ising antiferromagnets who referred to the excitations as a "Zeeman ladder" . These excitations have apparently been 
observed in CSC0CI3 |7j]. Such ladders of boundstates also occur in confining (1+1) dimensional quantum field theories 
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such as quantum electrodynamics (Q.E.D.) ||, the CP™ model || and a generalized "two-harmonic" version of the 
sine-Gordon model discussed below. Our approach to the spin-Peierls problem was partly inspired by the work of 
Coleman || on weakly coupled Q.E.D. 

We begin by considering a simple one-dimensional s=l/2 antiferromagnetic model without phonons: 



h = £[j(i + si-iy)^ ■ s l+1 + j 2 s t ■ s. l+1 }. 



(i) 



We first discuss the case 5 = 0. For J 2 < J 2c ~ .24 J the model is in a non-dimcrized gapless phase. On the other hand, 
for larger J 2 the groundstate is spontaneously dimerized and there is a gap jJJT^]. Thus, in principle, spontaneous 
dimerization could occur even without phonons. There is some evidence that the value of J 2 in CuGe03 may be close 
to J2C |0 but clearly phonons also play an important role in the spin-Peierls transition. We include J2 here for a 
different reason. In order to capture the essential physics of spontaneous dimerization without explicitly including 
phonons, it is neccessary to choose J 2 > J 2c - 

In the limit where J2 is only slightly larger than the critical value the low energy excitation spectrum can be 
determined by bosonization [fL2| . This gives the sine-Gordon model with Hamiltonian density: 



H = v 



1 



n- 



2 



dx 



g cos 2<p 



(2) 



with 8tt — 1 and g cx (J 2 /J 2c — !)■ Since this interaction has renormalization group scaling dimension (3 2 /4ir it is 
marginally relevant for J2 > J 2c leading to a soliton gap A s oc exp[— const/ (J2 — J 2c )]- It follows from the bosonization 
procedure that (f> is to be interpreted as an angular variable and that the discrete symmetry <f> — > (f> + n corresponds to 
translation by one site. Classically there are two groundstates at (f> = i 7r /2, corresponding to spontaneously broken 
translational symmetry. The s and s interpolate between these groundstates and have spin S z — ±1/2 depending 
on whether <f> rotates clockwise or counterclockwise with increasing x. These groundstates may be simplistically 
pictured as consisting of nearest neighbour dimers in one of the two possible patterns and the s or s as being a single 
unpaired spin separating the two different dimer patterns, as shown in Figure 1. The (presumably) exact results 
on the sine-Gordon model indicate that there are no other excitations besides the solitons and anti-solitons (and of 
course their multi-particle states) for this range of (3 fl3[| . In particular, unlike in a perturbative treatment in /3, there 
are no approximately harmonic excitations in addition to the topological ones. This is presumably a rather special 
feature of this particular model. The soliton width and spin correlation length obey £ s oc v/A s , diverging at J 2c . As 
we increase J2 further this soliton width decreases. At the special point J 2 — J/2 the exact groundstate is given by 
nearest neighbour dimers. In this entire region the translational symmetry is spontaneously broken so it follows from 
general principles that the excitation spectrum must contain solitons. 

— — * — — * — 

s s 

FIG. 1. An ss pair. The solid line represents a dimer singlet. 

Next we consider the model of Eq. (Q) with a small non-zero 5 orresponding to a mean field from the neighbouring 
chains preferring one of the two dimerization patterns. The staggered term leads to an additional sine-Gordon 
interaction: 



n 



ff 2 



n 2 



2 



dx 



g cos 2<f) + cS sin < 



(3) 



where c is a constant of 0(1). The spectrum is fundamentally different depending on the value of J 2 . For J 2 < J 2c we 
may ignore the marginally irrelevant cos 2<fi interaction in Eq. (|J) arising from the uniform terms in Eq. (0) . On the 
other hand, the sin0 interaction has dimension 1/2 and thus produces a gap oc S 2 ^ 3 (up to logarithmic corrections). 
The exact spectrum consists of a triplet and a higher energy singlet with A1/A3 = y/3 jl3]fj]. 

For J 2 > J 2c , there is a competition between spontaneous and explicit dimerization, represented by the g cos 2(f) 
and o5sin</> terms respectively in the sine-Gordon Hamiltonian. If 8 is very small then, while the true groundstate 
has cf> = —ir/2 the other state at d> = ir/2 has only slightly higher energy. Consider a ss configuration where the 
true groundstate (— n/2) occurs at x — * ±00 but the unstable groundstate occurs in between the ss pair, separated 
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by a distance x. Defining the s as the kink with <f> = —ir/2 at x — > — oo and the s as the kink with <f> = ir/2 at 
.x — » — oo, we see that the s must always be to the left of the s. If a; is much greater than the soliton width £ s then 
the classical energy of such a configuration is simply 2A S + c5x. This is true regardless of the spins of the two solitons 
(i.e. whether <f> winds clockwise or counter-clockwise). This linearly confining potential is crucial to the physics of 
spin-Peierls systems or indeed of any system with a competition between spontaneous and explicit breaking of a 
discrete symmetry. 

We now wish to treat this model quantum mechanically, in the small S limit. As far as we know, this "2-harmonic 
sine-Gordon model" model is not integrable. However, for small 5 we may follow Coleman's treatment of weakly 
coupled (l+l)-dimensional Q.E.D. p|. Consider a single ss pair. We expect it to give a spectrum of boundstates, 
due to the linear potential. Thus we write an effective Hamiltonian for the center of mass co-ordinate, x, of the pair: 

H = -M- S dx^ + C6x > (4) 

where M s w A s /v 2 is the soliton mass. We restrict x > and impose a vanishing boundary condition on the 
wavefunction at x = 0. The most important point about the eigenstates of this Hamiltonian is that they consist 
entirely of boundstates; there is no continuum. This follows because the potential keeps on increasing for all x so the 
soliton and anti-soliton never escape to infinity. A good idea of the nature of the spectrum can be obtained from the 
WKB approximation. The rfi 1 eigenvalue is given by: 



x 



dxyj M S {E B - c5x) = nw, (5) 



where xq is the classical turning point, xq = E B /cS. Thus the number of states with energy less than E B is: 



N(E ) « v sv > . (6) 

C061T 

We see that the density of states per unit energy diverges as 6 — ► 0. This is another peculiarity of a linear potential. 
As the strength of the potential goes to the free particle limit is obtained by the boundstates becoming more and 
more dense until they fill in the ss continuum. We have attached the superscript on E B to remind the reader that 
this is the binding energy of the ss pair; the total energy consists of this binding energy together with the ss rest 
mass energy: 

E n = 2A s + E B . (7) 

Actually, the results of the previous paragraph are only valid for low energy boundstates with E B < < A s where the 
non-relativistic approximation to the soliton dispersion relation may be used. Coleman extended the validity of this 
result by using the relativistic version of the WKB approximation. This is certainly valid for Q.E.D. We also expect it 
to be valid for our antiferromagnetic chain in the limit where c5 << A s << J, in which the theory is approximately 
Lorentz invariant up to the energy scale A s . Thus the free soliton energy may be written: 



E(p) = ^/A 2 +v 2 p 2 . (8) 

The WKB condition now becomes: 

J dpdx6{E n - 2E(p) - cSx) = i J dp[E n - 2E{p)]9[E n - 2E{p)] = 2nn. (9) 

Here E n is the full relativistic energy of the ss pair including both kinetic energy of the individual solitons and binding 
energy. So far this discussion ignores many body effects. In the weak coupling limit, as Coleman observed, these 
simply truncate this boundstate spectrum at E < 4A S . Any ss boundstatc of higher energy is unstable because it can 
decay into a pair of boundstates. If we imagine trying to pull an ss pair apart to oo it eventually becomes energetically 
favourable for "pair production" to occur so that we end up separating to oo two ss pairs. Setting E — 4:A S in Eq. 
(^) gives the number of stable boundstates: 

N .684A s /c<5. (10) 

We note that this 1 /S behavior is independent of the free soliton dispersion relation (although the prefactor depends 
on it). Shiba fq] encountered essentially a discrete lattice version of this Schroedinger equation in his work on Ising 
antiferromagnets, coming to similar conclusions about the spectrum. 
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So far we have ignored the soliton spin. The Hamiltonian of Eq. (0) is spin-independent so both the s and s can 
independently have spin up or down. Thus each of these boundstates corresponds to a degenerate triplet and singlet. 
When the s and s are close together (on the scale £ s ) their interaction will be a good deal more complicated; in 
particular it will be spin dependent. Fortunately, for very small S they stay far apart even in the lowest boundstate. 
Note that the n = 1 classical turning point is at xo oc 1/S 1 ^ 3 . As we increase S we expect the number of boundstates 
to decrease and the degeneracy beween triplet and singlet boundstates to be lifted. The set of stable boundstates 
must always lie below the two boundstate continuum in the coresponding spin channel. Clearly free solitons can never 
appear in the spectrum for non-zero 6. Eventually, when 5 >> (A s /J) 3 / 2 (but still S « 1) we expect to recover the 
spectrum of the J2 < .he model with one triplet and one singlet with the ratio A1/A3 = y3. It would be interesting 
to study the spectrum of the spin model numerically to test these ideas. (The expected behavior for J2 < J2 C and 
5 « 1 was recently confirmed H].) 

We now include phonons in our model, taking the Hamiltonian: 

H = J2l(J + auJSi ■ S l+1 + nf /M + [K/2)u* + si-iym]. (11) 

i 

Here Ui is the change in separation of two neighbouring ions from its uniform value and is the corresponding 
conjugate momentum; M is the ionic mass. (We have omitted acoustic phonons for simplicity. They presumably 
are not important for the spin-Peierls effect.) 8 represents a mean field arising from the coupling to neighbouring 
chains which favors one of the two possible lattice distortions. We may also keep a second nearest neighbour exchange 
coupling; it doesn't change the discussion qualitatively. While this Hamiltonian is considerably more complicated 
than the spin-only one considered above, we expect many features of the previous discussion to carry over. When 
<5 = we expect two degenerate groundstates with < Ui >= ±uq(— 1)\ The excitation spectrum will include s=l/2 
solitons and anti-solitons. Note that all excitations now involve both spin and phonon degrees of freedom. In addition 
to solitons, the excitation spectrum will presumably include other excitations corresponding to optical phonons. In 
principle there might also be integer spin magnetic quasi-particles below the two-soliton continuum. However, since 
the interaction between s and s vanishes at long distances for 6 = 0, no boundstates need occur in this limit, as in the 
previous spin-only model. Turning on a small 6 will again confine the solitons into ss pairs. In the small <5 limit we 
expect the number of stable boundstates to be proportional to 1/(5. Presumably spin ss boundstates can mix with 
optical phonons, which may reduce the number of stable s=0 boundstates. As 6 increases the number of boundstates 
will decrease. Free solitons can never appear in the spectrum. In addition to the stable boundstates there will also be 
two boundstate continua in the various spin channels (and also presumably a two phonon continuum). Evidently all 
boundstates must lie below the continuum in the corresponding spin channel in order to be stable. Numerical results 
on such a one-dimensional spin-phonon system would be highly desirable although considerably more difficult than 
for a spin-only system. 

Finally, let us consider a full three dimensional spin-Peierls Hamiltonian. A simple model would consist of chains 
with the Hamiltonian of Eq. ([ll]) (with 5 — 0) together with an inter-chain phonon coupling: 

^'EE u iR u *R-> ( 12 ) 

' <R,R/> 

where R labels the lattice points perpendicular to the chains. In the ordered phase < u.^ >= (— l) l ito- Treating 
the interchain coupling in mean field theory gives the one dimensional model of Eq. ( [Tif ) with S = ZK'uq where 
Z is the chain co-ordination number, uq should then be determined self-consistently by solving the one-dimensional 
model. This is essentially the approach advocated by Khomskii et al. H. Note that it is a rather standard approach 



to various quasi-one-dimensional systems 15 . In the disordered phase, 8 = so solitons can exist as independent 
excitations. In the ordered phase 5 > and free solitons cannot occur due to the confining potential. While based on 
a mean field treatment of inter-chain couplings this conclusion is presumably much more general. It essentially follows 
from Landau's argument that any non-zero density of free solitons leads to the destruction of long range order in a 
one dimensional system. The theoretical discussion in Ref. || essentially considered a purely one-dimensional spin- 
phonon model which does indeed contain free solitons. The magnon was regarded as an ss boundstate that could 
occur in this model. However, it is not permissible to ignore inter-chain elastic coupling in the dimerized phase. This 
coupling at the same time stabilizes the dimerized phase up to a finite crticial temperature and confines the solitons. 
Note that in the small K' limit, the energy scale of all excitations is set by the one-dimensional model and therefore 
the magnon gap should be given by twice the soliton gap (assuming no boundstates in the one-dimensional model for 
(5 = 0). Even if K' is not very small, its effects become smaller near a transition into a non-dimerized phase, driven 
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by temperature, field or impurity concentration. Thus, at least naively, one might expect the present approach to 
become more valid near such a transition. 

What does this approach tell us about CuGeC>3? We may interpret the observed magnon as a spin 1 ss boundstate. 
Likewise, if there is a stable singlet quasi-particle, as suggested by Raman scattering Jl6| , we may interpret it as an 
s=0 ss boundstate. As remarked above, there appears to be no sharp distinction, in general, between such an s=0 ss 
boundstate and an optical phonon in a spin Peierls system. In Raman scattering the photon couples both to lattice 
displacements and to Si ■ Si+i- Such spin excitations might also be observable in neutron scattering because they 
couple to lattice displacements which are excited by neutron scattering from the ionic nuclei. One precise conclusion 
from the present approach is that there can be no soliton continuum in the spin-Peierls phase. The continuum in the 
magnetic neutron scattering cross-section neccesarily starts at twice the gap to the lowest spin triplet ss boundstate 
corresponding to a 2-magnon continuum. This is consistent with the data presented in fl. The apparent absence of 
additional ss boundstates in CuGe03 (besides the magnon and possibly one singlet) can presumably be attributed 
to the relatively large value of the inter-chain phonon coupling, K' . Likewise the failure to observe an ss continuum 
above T$p can be so attributed since then T$p is of the order or greater than A s , smearing the continuum threshold. 
It would be interesting to find spin-Peierls materials which were more highly one-dimensional with respect not only 
to their magnetic exchange couplings but also their elastic couplings. In such materials additional ss boundstates 
should exist below T$p and the ss continuum could perhaps be observed above T$p. One way in which enhanced 
one-dimensionality might occur is in a system where the next nearest neighbour Heisenberg coupling J2 > J2 C ~ .24 J. 
In this case the soliton gap could be determined primarily by the magnetic exchange energies and might be large 
compared to K' or other phonon energy scales. 
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